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We present a next-to-leading order QCD calculation for single-inclusive high-pj- jet pro- 
duction in longitudinally polarized pp collisions within the "small-cone" approximation. 
The fully analytical expressions obtained for the underlying partonic hard-scattering cross 
sections greatly facilitate the analysis of upcoming BNL-RHIC data on the double-spin 
asymmetry for this process in terms of the unknown polarization of gluons in the 
nucleon. We simultaneously rederive the corresponding QCD corrections to unpolarized 
scattering and confirm the results existing in the literature. We also numerically compare 
to results obtained with Monte-Carlo methods and assess the range of validity of the 
"small-cone" approximation for the kinematics relevant at BNL-RHIC. 



Abstract 



I Introduction 



The first successful runs of the BNL-RHIC with polarized proton beams mark a new era in 
spin physics. Very inelastic pp collisions at high energies open up unique possibilities to answer 
many interesting questions, first and foremost perhaps those concerning the still unknown gluon 
polarization in the nucleon, Ag p. In the near-term this prime goal can be accomplished by 
studying double spin asymmetries All for the inclusive production of hadrons and jets. Both 
are copiously produced at high energies, and luminosity requirements are rather modest for 
a first determination of Ag. Measurements of the unpolarized neutral pion cross section at a 
center-of-mass (c.m.s.) energy of = 200 GeV, for both central and forward pseudo-rapidities 
1] by Phenix [2^ and Star respectively, have shown good agreement with perturbative QCD 
(pQCD) calculations, even down to unexpectedly small values of the pion transverse momentum 
Pt of about 1.5 GeV. This boosts confidence that similar measurements with polarized protons, 
too, may be interpreted in terms of pQCD. Very recently, Phenix has published first results 
for All in PP vr'^X at moderate Pt, though preliminary and with rather limited statistics 
Surprisingly, the data show a trend towards a negative and sizable All which, if it persists, 
would appear to be impossible to accommodate in the framework of leading twist pQCD 
More data are clearly needed before any definitive conclusions can be drawn. It is expected that 
Star will soon publish data on the spin asymmetry Aj^^ the closely related jet production. 
It will be particularly interesting to see whether these data will show a similar trend as the 
Phenix tt^ data. This paper presents new theoretical pQCD calculations and predictions for the 
spin asymmetry in jet production, which appears timely in view of the experimental situation. 

In order to make reliable quantitative predictions and to analyze upcoming data in terms 
of spin-dependent parton densities it is crucial to know the next-to-leading order (NLO) QCD 
corrections to the lowest order (LO) Born approximation for the partonic hard-scatterings. 
NLO corrections can be sizable and may affect the spin asymmetries. More importantly, the 
dependence on unphysical factorization and renormalization scales is expected to be reduced by 
the inclusion of higher order corrections. In case of jet production, NLO corrections are also of 
particular importance since it is only at NLO that the QCD structure of the jet starts to play 
a role in the theoretical description, providing for the first time the possibility and necessity to 
match the experimental conditions imposed to group final-state partons into a jet. 

The past few years have seen much progress in calculations of NLO QCD corrections to spin- 
dependent processes relevant for the RHIC spin program. In particular for inclusive hadron 
production, NLO results have been obtained recently both at an analytical level as well 
as using Monte-Carlo integration techniques [Zj. While the latter basically allow to compute 
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any infrared-safe cross section numerically, the former is restricted to inclusive hadron spectra. 
However, fully analytical calculations lead to much faster and more efficient computer codes 
as all singularities arising in intermediate steps have explicitly canceled and are not subject to 
delicate numerical treatments. Such calculations will greatly facilitate a future "global analysis" 
of RHIC data in terms of polarized parton densities, e.g., along the lines described in jH]. 

The aim of this paper is to obtain (approximate) analytical NLO QCD results also in 
the case of single-inclusive high-p^ jet production in longitudinally polarized pp collisions. 
NLO corrections to polarized jet production have been determined before, using Monte-Carlo 
methods 9J. However, for the reason just given, we believe that our calculation is a very 
useful addition and of great relevance in practice. Its results will be immediately usable in a 
phenomenological analysis of forthcoming Star data. 

The actual calculation of the partonic cross sections up to NLO is a formidable task. Fortu- 
nately, inclusive hadron and jet production proceed through the same set of partonic subpro- 
cesses, and it turns out that we can make use of much of what we already calculated previously 
in the case of hadron production 6 . The main difference is in the treatment of final-state 
singularities: in case of inclusive hadron production the phase space of all unobserved partons 
is fully integrated over. This leads to final-state singularities which are to be absorbed into the 
parton-to-hadron fragmentation functions. However, a jet is essentially defined as a transverse- 
energy deposition within a certain cone centered around pseudo-rapidity rf'^^ and azimuth 
For such an observable, all final-state singularities cancel ^Ul- As we will see below, despite 
this difference, it is possible to transform the single-parton-inclusive cross sections that we have 
calculated earlier jH] into single-inclusive jet cross sections. This can even be done at an an- 
alytical level if one assumes that the jet cone is rather narrow ("small-cone approximation", 
SCA). This is the strategy we will use in this paper. We note that the SCA was first introduced 
many years ago fOlE] and has been applied in computations of unpolarized single-inclusive jet 
cross sections ^2]- We have re-done the calculations in the unpolarized case and fully confirm 
the results available in the literature ^2]- By comparing our polarized cross sections obtained 
within the SCA to results of the NLO Monte-Carlo "parton generator" of we can assess 
the range of applicability of the SCA. It turns out that for the kinematics relevant at RHIC 
the SCA is very accurate, and that our numerical code is orders of magnitude faster than the 
Monte-Carlo one. Thus our results will indeed be useful for the analysis of upcoming data in 
terms of polarized parton densities. 

The outline of the paper is as follows: after defining our notation, the next section is mainly 
devoted to present in some detail the necessary technical framework to convert single-parton- 
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inclusive cross sections into jet cross sections. In Sec. Ill we present some phenomenological 
applications of our results, all tailored to upcoming measurements with the Star experiment 
at the BNL-RHIC. First we compare the SCA to the Monte-Carlo results obtained for various 
cone sizes. Next we examine the size of the NLO corrections, the reduction of dependence on 
unphysical scales in NLO, and the sensitivity of the double-spin asymmetry All to ^g. We 
summarize the main results in Sec. IV. 

II Jet Production in the Small-Cone Approximation 
II. 1 Jet Definition and SCA 

We will consider the single-inclusive jet cross section for the reaction pp — > jet X which counts 
all events with a jet of a given transverse energy E^^^ and pseudo-rapidity r^-''^*. Cross sections are 
infinite unless a finite jet "size" is imposed. Jets are not intrinsically well defined objects, but 
may be constructed in somewhat different fashions from a set of close-by final-state particles. 
In this work we define the jet four-momentum as the sum of the four-momenta of all particles 
inside a geometrical cone of half-aperture 5 around the jet axis, given by the three-momentum 
of the jet, see Fig. [T] We treat the jet in the SCA jTU [T2] as this allows for a largely analytical 
calculation of the cross section. 

Other jet definitions and algorithms are clearly possible and, in fact, more widely used in 
practice [13-17]. A popular choice P^, also adopted in the jet analyses at STAR, is to define 
a jet as a deposition of total transverse energy E'^ of all final-state particles that fulfill 

(r/j^* - T]'f + - <p'f < . (1) 

Here 77* and 0* denote the pseudo-rapidities and azimuthal angles of the particles, and R is the 
jet cone aperture. The jet variables are defined by i?;^* = Ylii^Ti V*^* — — ln[tan(^-''^Y2)] = 
J2iE'TVyE}^\ and = ^.E*.0VE^^\ It has been shown fW^ that in the SCA the jet 
definition we have adopted becomes equivalent to that given by Eq. ((T)), provided 6 is chosen 
as R/ cosh (77^*^*) . We have verified this equivalence numerically using the Monte-Carlo code of 
jU]. Differences arise only at and are negligible for small cone apertures. 

The SCA may be viewed as an expansion of the partonic cross sections around 5 = (or, 
equivalently, R = 0). For small 6, the dependence on the cone size is of the form Alog{6) + B + 
In this work we determine the coefficients A and B at an analytical level. We neglect 
the 0{S'^) pieces; we will demonstrate in Sec. Ill that this is a surprisingly good approximation 
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Figure 1: Sketch of single-inclusive jet production in p{Pa)p{Pb) ^ jet(Pjct) 

even for (experimentally relevant) cone sizes of up to i? ^ 0.7. This observation was also made 
some time ago for unpolarized jet cross sections by comparing the SCA to a calculation where 
the SCA was extended numerically to finite cone sizes [T^ l2Uj . 



II. 2 Notation and Outline of the Calculation 



According to the factorization theorem for high-py processes j2T], the spin-dependent cross 
section for single-inclusive jet production p{Pa)p{Pb) jet(Pjet) X can be written as 
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where the superscripts in Eq. (jSj) denote the helicities of the colliding polarized protons. The 
sum in Eq. (jS)) is over all contributing partonic channels a -|- 6 — > jet -|- X, with their associated 
LO and NLO partonic cross sections dAa^J^-^^^ and dAa'^jj_^^^^^^, respectively. The latter are 
defined in complete analogy with Eq. Q, the helicities now referring to partonic ones. In 
Eq. (jHI) we have introduced the dimensionless variables V and W which are defined in terms of 



pi^^ and T]^^^ as 
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jet 

Ft 
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with S = [Pa + PbY the available hadronic c.m.s. energy squared. The corresponding parton- 
level variables read 

t —u 

V = l + - , W = — , S = {j)a+ Vhf , t={pa- Pjetf , U = {pb - Pjctf , (5) 

S S + t 

where Pjet is the four-momentum of the jet {P^et — in the SCA). The information on the spin 
structure of the proton in Eq. (jH} is contained in the spin-dependent parton densities A/a ^. 
They are probed at momentum fractions Xa,b given by 

..^^ and ..^i^ . (6) 

VW 1 — V 

The factorized structure of Eq. (jH)) dictates the appearance of the factorization scale fip which 
is of the order of the hard scale in the reaction, p^^*, but not further specified. The same 
is true for the scale {Ir associated with the renormalization of the running strong coupling 
c^sifi'R)- Compared to single-inclusive hadron production E] the jet observables in Eq. Q 
have the advantage of being free of uncertainties associated with the factorization of final-state 
singularities into non-perturbative parton-to-hadron fragmentation functions at a scale fi'p. 
Hence, jet production at RHIC is expected to provide particularly clean and useful information 
on the spin structure of the proton. 

As indicated above, the jet cross sections at the parton-level may be evaluated in pQCD as 
an expansion in the strong coupling a<j. 



ab->jctX + ^ "''-^^ab- 

The LO approximation of Eq. ((Tj), rfAo-^^^^j^^j,^, is obtained from evaluating all basic 2 — > 2 QCD 



dMab-.-,.tx = dAaZli.^^ + ^c^AaiiLjetx + • • ■ • (7) 

.(0) 



scattering diagrams, that is, X consists of only one parton recoiling from the other parton 
producing the observed jet. Contrary to single-inclusive hadron production where different 
final-state partons have to be weighted with different fragmentation functions leading to ten 
separate LO channels [H], for a jet cross section no distinction is made between different quark 
flavors or gluons producing the jet. All processes with the same initial-state partons have to 
be summed appropriately, which is precisely the reason why final-state singularities cancel at 
higher orders. One then has only six different subprocesses : 
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Figure 2: Contributions to the single-inclusive jet cross section from partonic reactions where only 
one parton is in the cone (left hand side) and where two essentially collinear partons, j and k, form 
a narrow jet (right hand side). 

For the computation of the NLO corrections to Eq. ((7j), dAa^^ij^-^^j^, we would hke to make 
use as much as possible of the single-parton-inclusive cross sections that we have calculated 
earlier [0] . In the following we will demonstrate how to convert analytically the single-parton- 
inclusive cross sections into jet cross sections. In order to achieve this, we use the SCA and 
that the jet is formed either by a single final-state parton or by two partons essentially collinear 
to each other, see Fig. |21 

The connection between the two types of inclusive partonic cross sections can be best 
understood in the following way fn\ 120] : let us consider the generic 2^3 partonic reaction 
0'{Pa)b{pb) j{Pj) k{pk) l{pi)- Partons a and b are longitudinally polarized here, which is 
however not really of any relevance in what follows. In we have already computed at an 
analytical level the single-parton inclusive cross sections for ab jX, ab kX, and ab —>■ IX. 
Each one of these already contains a full sum over unobserved partonic final states X, that is, 
ab jX would consist of [ab —>■ j{kili)] + [ab —>■ j(/i;2^2)] + • • •, in case several different final 
states are allowed. For instance, the channel qg - 
reactions qg -> q{gg), qg q{qq), qg g(g'g')' 
parton processes, ab — > jX, ab — ^ kX, and ab 
production, for which they have to be added. 



■ qX receives contributions from the partonic 
With these stipulations, all three inclusive- 
IX, contribute also to single-inclusive jet 



Looking closer at the cross section for ab — ^ jX one has to distinguish configurations where 
only parton j is in the cone and forms the jet from situations where, for example, partons j 
and k (and similarly, j and /) are both in the cone, see Fig. El In the latter case, the single- 
parton-inclusive cross section still assumes that only parton j is observed, whereas for the jet 
cross section we want to define the jet from both partons inside the cone. We therefore subtract 
those parts from the single-parton-inclusive cross section for which j is observed, but k (or /) 
are in the cone. We proceed, of course, in the same way for the cross sections with k or I being 
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the observed parton. Finally, we have to add pieces for which j and k (or j and I, or k and 
I) are both in the cone, and for which j and k are forming the jet together. Note that these 
contributions are symmetric under exchange of j and k. 

To be more specific, let us denote the longitudinally polarized single-parton-inclusive cross 
section for ab jX by dAaj, the one where still j is observed, but k is also in the cone by 
(iAo"j(fc), and the one where j and k are in the cone and form the jet together by dAajk- We 
then have the following final expression for the desired partonic jet cross section: 

dAaab^jetx = [dAcTj - dAaj(k) - dAaj(i)] + [dA&k - dA&k^j) - dA&k^i)] 

+ [dAai - dAai(j) - dAa^k)] + dAajk + dAaji + dA&ki ■ (9) 

All final-state singularities of the individual cross sections are guaranteed to cancel in Eq. (jU)). 
Technically, because the dA&j have already been made finite by an MS subtraction ^ , we will 
first compute the combination 

- dAaj(^k) - d^^jH) - dA&k^j) - dAdkii) - dA&K^j) - dAanj,) + dAdjk + dAdji + dA&ki , (10) 

which will have final-state coUinear singularities, the same (but with opposite sign) that were 
originally present in dAaj + dAdk + dAdi. To obtain a finite answer for the combination an 
MS subtraction then has to be performed. A powerful check for the correct implementation of 
all contributions is to see that the dependence on the final-state factorization scale ^'p present 
in the intermediate results disappears identically for each of the subprocesses listed in Eq. (jHJ 
when dAdj + dA&k + dA&i, taken from is finally combined with (fTIH) . 

Clearly, the pieces dAdji^k) and dAajk that we will subtract and add are dominated by con- 
figurations with two collinear final-state partons. These configurations produce the A\og{6) + 
B + (9(5^) behavior in 6 mentioned earlier. For collinear kinematics, the calculation dramat- 
ically simplifies because the 2 — > 3 matrix elements then factorize into 2 —>■ 2 ones and LO 
splitting functions, allowing the calculation to be done largely analytically. The narrower one 
chooses the cone, the more dominant will the "collinear" terms become. This is the reasoning 
behind the use of the SCA. 

We note that the separation in Eq. © in the context of the SCA was already considered in 
Refs. jTTl 1201 unpolarized case. In these references also some details of the calculations 

of the dcTji^k) and ddjk were given. The aim of the following subsections is to make this paper 
self-contained by providing all essential calculational details, without however going to excessive 
length. In the next two subsections, we will address dAaj(j^) and dAajk separately. 
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II. 3 Calculation of dAaj(k) 



In this section we show how to calculate the "subtraction" contributions dAaj^k) to Eq. (0) 
when parton j alone forms the jet but parton k is also in the cone. We note that this part of 
the calculation is very similar to what was done for the isolated prompt photon cross section 
[22]; except for the fact that all possible parton-parton splittings will occur in our case. The 
leading contributions in the SCA result from a branching of an intermediate parton I (which 
could, of course, be of the same type as j and/or k), see Fig. |21 In the coUinear approximation, 
the spin-dependent matrix element squared A\M\li^^ji^i for the process ab jkl factorizes into 



1 



2pj • Pk 



vw 



vw 



PJ}{z = 1 - V + vw) , (11) 



where all quantities are in d = 4 — 2e dimensions in order to regularize coUinear singularities. 
The superscript "<" denotes that the d-dimensional I ^ j unpolarized splitting function Pji{z) 
is strictly at 2; < 1, that is, without its 5{1 — z) contribution if / = j. A|M|^^_„ denotes the 



spin-dependent matrix element squared for the 2 

2pj ■ Pk 



9 V il — w) 



2 process ah 
■(1 - cosOjk) 



II. In addition, we find 



(12) 



1 — V + vw 

with 9jk the angle between partons j and k. The factor l/(2pj -pk) in (fTTj) constitutes the only 
dependence on 9jk in the coUinear approximation. Integrating it over the 2 — 3 phase space 
gives 
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(13) 



where 6 is the cone opening, v' is as in Eq. (fTT|) . and the factor in square brackets in the first 
line is the usual 2^2 phase space in d dimensions. For small 6 the angular integral in Eq. (fT^ 
is readily evaluated as 
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(14) 



For the d-dimensional partonic cross section at 0{a^g) in the collinear approximation we then 
arrive at 

dAdab^n 
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dvdw 
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Pjj{z = 1 — f + vw) 



Til-e)l-v + vw 
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2n 



(15) 



Here, dAaab^n/dv denotes the ci- dimensional Born cross section for ab —>■ II. One can observe 
in Eq. (fTH|) the expected separation into a universal part associated with the splitting function 
Pj7, and a part related to an underlying 2^2 partonic scattering. 

From now on we will mainly focus on the splitting part, keeping in mind, however, that the 
appropriate sum over all partonic channels and splittings is implicitly understood in the end. 
The further evaluation is particularly simple when the splitting happens to be non-diagonal, 
i.e. / 7^ j. In that case, the full splitting function Pji is regular at 2; = 1, i.e., at w = 1, 
and all terms may be safely expanded in e. The collinear singularity arising in Eq. (|15p was 
also present in the single-parton-inclusive cross sections dA&j, where it was subtracted at a 
scale ^'p into the bare parton-to-hadron fragmentation functions in the MS scheme |6j. This 
factorization is not appropriate for a jet cross section. In order to compensate for it, we have 
to apply the same subtraction also to Eq. (jisp . Any dependence on the factorization scale fi'p 
will then drop out for the inclusive jet cross section in the sum ^ over partonic scatterings, 
as discussed above. The appropriate factorization term for Eq. (fT3jl reads 

C?A(Tfact as dAaab-,11 ( I vw \ (4) 

— ' V = ^ji = 1 — f + vw) 



dvdw 271 dv \ 1 — f + vw J 

. (16) 



e J 1 — V + vw 

where now Pjf is the usual four-dimensional splitting function. Expressing the d-dimensional 

one as 

P,r{z) = P^t\^)+ePij\z) , (17) 
and adding Eqs. (fTKj) and (PT)|) we arrive at the final result for dA&j^k) if j 7^ I, 

dAd-j(k) _ ots V dAa[lz^]i [ , _ vw 



V 



dvdw 27r 1 — f + vw dv \ 1 — v + vw 

p')f{z =1-v + vw)\yi -p^ -J -Pfj\z = l-v + vw) 

which is finite and shows the expected logarithmic dependence on 5. 

For a diagonal splitting function in Eq. (jl5|) . i.e., j = I, the situation is somewhat more 
complicated due to the infrared singularity of Pjj{z) at 2; = 1. The singularity is regularized 
by the factor (1 — w)^^'^ in (fT3j) and gives rise to a l/e"^ pole. This pole is then canceled by 
contributions from dAajk, i.e., the cross section with both partons in the cone forming the jet. 
Rather than giving the lengthy (regularized) expression for the diagonal case, we therefore first 
turn to the calculation of dAajk- 
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II. 4 Calculation of dAajk 



Again, the leading configurations are those for which an outgoing intermediate parton / splits 
collinearly into two. We now have to deal with partons j and k in the cone producing the jet 
with four-momentum Pjet = Pj + Pk- Let us first note that for the dimensional 2 — > 3 phase 
space one finds 



(27r)5-4£ sv ^ 7 ^ Ek 

The cross sections dAajk thus contribute only at w = 1, i.e., with 2 — > 2 kinematics. As 
discussed in the previous subsection the only dependence on 6jk, the angle between the two 
partons j and k, stems from the propagator oc l/{2pj -pk)- Integrating this term over the phase 
space (fTn|) one finds after some algebra: 



dPS, 1 - . - 



1 f^nV 1 . 

v[l — v)\ 



TT \ s J T{l-e) 



7? /F^X-"^ /-emax „• l-2£/l 



dvdw (2pj ■ Pk) 

r{l-e)Jo ' El\s ) Jo n-cos%fc ^ ^ 
where 6j is the angle of parton j with respect to the direction of the jet, Pjet- 

In order to find the relation between 6j and 6jk we first write 

n _ Pj ■ Pk _ • (^j^^ ~ _ li^jctl COsdj - Ej 

" b-.l \Pk\ ~ E, (Eje, - E,) - Ejet - E, 

and then use 

^jet = (Pj + Pk? = E] + (^jet - Ejf + 2^,(^jet - Ej) COS 9jk • (22) 

Solving Eqs. 1)211) and ()22|) for ^j, one finds in the coUinear limit 

e, ^ ^^^^e,k . (23) 

In the same way, one finds for the angle 6k between parton k and the direction of the jet 

Ok ^ 4^e^k . (24) 

It is important to keep in mind that neither of the two particles j and k is allowed to be outside 
the cone. In other words, the upper limit 6'max in Eq- fp^ in terms of the cone opening 5 is 

Ek > Ej : 6'max = S , 

Ej > Ek : ^^ax = . (25) 

Ej 
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We note that the integration over 9j in Eq. (j'iUp could also be written as an integration over 
the invariant mass \JP^t of ^^e jet. Consistently with the SCA, we have previously neglected 
the jet mass everywhere. However, the singularity of the cross section at Ojk = 0, where the jet 
mass vanishes, requires us to integrate over a finite region < Pj^^ < 0{S'^E?^^) here. 

We finally obtain, after performing the 6j integration in Eq. ()20|) and introducing the variable 



dPS; 



dvdw {2pj ■ pk) 



1 UnV 1 
8^ VTy T{l-e] 



Hi-v)]- 



87r2 



6(1 -w) 



X : 



r 1 



d^ 



0"'^ 0(^-2 



(26) 



Here, 6 denotes the Heaviside step-function, delineating the two regions specified in Eq. (j^SJ. 

Eq. (Illj) obviously continues to hold. The difference between the calculation of dAaj(^k) in 
Sec. II. 3 and of dAcjjk here is that we now need to integrate over the argument of the splitting 
function appearing in Eq. ()1H) . as ^ in Eq. (j^UI) plays the role of the momentum fraction z. To 
do so, let us denote 



di 



Then one straightforwardly computes the four relevant integrals J„ 
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(27) 



^99 



1 

2 

2C, 



+ e I -I + 



TT 



31n(2) 



'■91 ' 



1 11 

e 6 



137 _ n 

36 ^ 3 3 ^ 



:28i 



with Ca = 3 and Cp = 4/3 the usual SU(3) Casimir operators. Note that the 1/e pole terms 
in Eqs. fl28|) . along with the overall 1/e singularity in Eq. (f^ . give double poles that exactly 
cancel the ones arising for collinear sphttings in the dA&j(^k) (see the discussion at the end of 
Sec. II. 3). This cancelation will be demonstrated in the next section. Furthermore, the term 
—3/2 in the first line of Eqs. (j28|) is identical to the negative of the coefficient of the 6{1 — z) in 
the full Pqq splitting function. It is extremely important that we obtain this term: clearly, the 
factorization subtraction in Eq. (fTBj) contains the full splitting function, whereas all collinear 
limits of 2 ^ 3 matrix elements, Eq. (fTT|) . give only the splitting functions without their 6{l — z) 
contributions. Without the term —3/2 above we would not be able to cancel the final-state 
singularities. The same happens if the intermediate parton / is a gluon: an outgoing gluon 
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may split into a gg or a qq pair, the latter coming in Nf "active" flavors. Again we flnd in 
Eqs. ()28|1 the necessary —11/6 and 2/3 factors, which combine to the LO QCD /5-function, 
Po = llC^/3 — 2Nf/3, appearing in the 5(1 — z) contribution to Pgg. 



II. 5 Collecting Results and Cancelation of Singularities 



The rest of the calculation is essentially an exercise in bookkeeping of all of the many interme- 
diate terms appearing in the various contributions to Eq. Q. 

To be more speciflc, let us sketch an example: a process with a qqg final state. We already 
know dAcTq, dA&q, and dAdg from the inclusive hadron calculation [H]. The pieces dAag(q) and 
dAag(^q) rely on the splitting function Pgg and are treated by Eq. (fTK|) . Likewise, dAdq^q) is 
related to Pgg. One only has to attach the appropriate LO cross section [depending on v' as 
defined above in Eq. (fTT|) ]. On the other hand, dAag(^g) and dAdqi^g) depend on the diagonal 
splitting function Pgg and are singular. It is therefore convenient to combine them directly with 
dAdqg and dAdgg, respectively. The combination dAaqi^g) — dAaqg is, after MS subtraction of 
final-state coUinear singularities at a scale fj,'p, proportional to 



dAdqi^g) — dAa, 
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(29) 



This result is completely finite and may serve as a "building block" whenever there is an 
outgoing quark radiating a gluon. Only the appropriate longitudinally polarized LO cross 
section dAdah^qil dv has to be attached to Eq. (|^. summed over all final-state partons /. The 
"+" -distribution in Eq. is defined as usual by its integration with an appropriate test 
function: 



dw [g{w)] f{w) 



dwg{w) [/(«;) -/(I)] . 



(30) 



The final piece in this example is dAoqq which involves the integral Iqg in Eq. (j28jl and 
therefore is also singular. It is crucial here that we sum over all possible final state channels for 
two given colliding partons: if there is a contribution involving dA&qq, there must also be one 
involving dAcjgg. With the appropriate combinatorial prefactors these two will combine to a 
/3o, as discussed above. Including also (iA(9"g(g), and after MS subtraction of final-state coUinear 
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singularities, one obtains for the "building block" for intermediate gluons, i.e. I = g, 



2dAa 



dig) 



2Nf dAa. 



qq 
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[l + (1 - + vw{2 -2v + vw)] In 



ln(l 



w 



w 



w 



I 2 

^^F 



(31) 



again finite. Recall that contributions from dAdqi^q), which also appear in Eq. 0, are associated 
with the non-diagonal splitting function Pqg and thus are finite by themselves and hence not 
included in (1311) . They are straightforwardly treated according to Eq. (fTHj) . Our main equations, 
Eqs. (|TH|) . (j2ni), and (jSH), all show the expected "logarithmic plus constant", i.e., ^log(5) + i3 + 
0((5^), dependence on the jet cone opening 5. We note that the terms oc [ln(l — w)/{l — 
in Eqs. fl29|) and (jHTj) . which are the leading contributions at — > 1, cancel identical pieces 
arising from the observed final-state parton in the single-parton inclusive cross sections dAdq, 
dAdg, so that such distributions are absent in the jet cross section. This finding is in accordance 
with results found in a large-w "threshold" resummation calculation of the jet cross section, 
when the jet is allowed to be massive at partonic threshold 

With these prerequisites at hand one can compute all relevant NLO partonic single-inclusive 
jet cross sections listed in Eq. (jSI) in the SCA. The final analytical results are too lengthy to 
be given here but can be found in a Fortran code which is available upon request from the 
authors. As mentioned already, a powerful check for the correctness of the calculation is the 
cancelation of any dependence on the final state factorization scale n'p when, according to 
Eq. Q, the single-parton- inclusive cross sections from jH] are combined with the appropriate 
combinations ofdAajf^k) and dAajk according to Eqs. (fTHj). (j^ . and (jSH). All of our final results 
pass, of course, this consistency check. Another important check of the procedure outlined here 
comes from the computation of the unpolarized jet cross section in the SCA - the building 
blocks in Eqs. and have no memory of the polarization of the colliding partons which 
only enters through the LO 2^2 cross sections attached to them. We fully agree at an 
analytical level with the results of which can be retrieved from their Fortran code, after 
an appropriate transformation to the MS factorization scheme. 
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Figure 3: Ratio of the NLO (a) unpolarized and (b) polarized inclusive jet cross section in the SCA 
and within the full Monte-Carlo approach of [9J for = 200 GeV, several bins in p}^^, and three 
different cone sizes R. 



Ill Numerical Results 



We now turn to a phenomenological study of single-inclusive high-pr jet production in polarized 
pp collisions at RHIC. Regarding our results obtained in the SCA, the most important question 
is, of course, how accurate the approximation actually is in cases of practical relevance. We 
will investigate this first, by comparing our results to those of the Monte-Carlo code of |9]. For 
this purpose, we choose the recent CTEQ6M set of unpolarized parton densities in the 
calculation of the unpolarized cross section, and the NLO "standard" set of GRSV 25^ for the 
polarized case. We assume the kinematical coverage of the Star experiment, which can detect 
jets in the pseudo-rapidity range —1 < rf^^ < 1, over which we integrate. 

In Fig. ini we compare the results of the full Monte-Carlo NLO jet calculation jH] to the 
results within the SCA for = 200 GeV, several bins in pj,^, and three different cone sizes R. 
The histograms on the left show the unpolarized results, the right ones give the comparison for 
the polarized case. It turns out that even for the rather large cone radius of R = 0.7 the SCA 
gives still very acceptable results within ten percent or less of the full Monte-Carlo calculation. 
For the unpolarized jet cross section such an observation was already made in ^niEDj- In the 
polarized case, however, the success of the SCA up to large R is non-trivial due to possible 
cancelations between the two helicity configurations in Eq. Not unexpectedly, for very 
large R = 1.0 the SCA starts to break down because neglected contributions proportional to 
0{R^) become relevant then. It is expected though that a cone size between 0.4 and 0.7 will 
be chosen by the Star collaboration in their forthcoming analysis; larger sizes are not really 
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Figure 4: Unpolarized and polarized inclusive jet production cross sections in NLO (solid) and LO 
(dashed lines) at (a) = 200 GeV and (b) = 500 GeV for cone sizes R = 0.4 and R = 0.7, 
respectively. In each case the lower panel shows the ratios of the NLO and LO results. 



practical in view of tlie limited angular acceptance of the detector. 

Fig. 121 demonstrates that our results based on the SCA are sufficiently accurate to be used 
in analyses of forthcoming data on jet cross sections and spin asymmetries from RHIC. We 
emphasize that numerically stable results for the full pip' spectrum can be obtained with our 
computer code very fast and efficiently, in a matter of minutes. This malces our code an ideal 
ingredient for future "global" analyses of RHIC jet data in terms of polarized parton densities. 
This is a clear advantage over a Monte-Carlo code with its huge numerical complexity, which 
yields results with rather large numerical fluctuations (still visible in the histograms shown 
in Fig. inj even after hours of running. This is even more true for the polarized cross section 
due to large cancelations between the two helicity configurations in Eq. (0). That said, our 
present calculation can only be used to describe single-inclusive jet cross sections, whereas the 
Monte-Carlo code is much more flexible concerning the observables that can be predicted. 

Heartened by the good agreement between our code and the full Monte-Carlo calculation, 
we will now present a few predictions for RHIC. We will be very brief here because many 
phenomenological results for jet production have already been presented in [D]. We focus on 
the most interesting questions: the importance of the NLO corrections, the residual dependence 
on the unphysical scales and in Eq. Q at NLO, and the sensitivity of the double spin 
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Figure 5: As in Figs. 111(a) and (b) but now showing the scale dependence of the polarized cross 
section in LO and NLO in the range p^^V^ < A^f = A^i? < '^Pr^- We have rescaled the LO results 
by 0.1 to separate them better from the NLO ones. In each case the solid lines correspond to the 
choice where both scales are set to p^^*. 



asymmetry 



dAa 
da 



(32) 



to the still unknown gluon polarization Ag. Predictions for in immediate demand for 

an extraction of Ag at RHIC in the very near future. In comparison to P, there are also some 
new results in our analysis. 

Figure |3] shows our results for the unpolarized and polarized p^^*-spectra of single-inclusive 
jets at NLO and LO, integrated over -1 < r/J^* < 1, for = 200 GeV and = 500 GeV and 
cone sizes R = OA and R = 0.7, respectively. We have set the scales ^r = = Pt"- Again we 
have used CTEQ6M parton densities |^ in the unpolarized case and the GRSV "standard" 
set for the polarized cross section, always performing the NLO (LO) calculations of cross 
sections using NLO (LO) sets of parton distribution functions and the two-loop (one-loop) 
expression for ag- The lower part of the figure displays in each case the so-called "i^-factor" 

c/(A)aNLO 



K 



(33) 



d{A)a^o ■ 

One can see that the NLO corrections are fairly moderate and of similar importance in the 
polarized and the unpolarized cases. 

Figure shows the scale dependence of the spin-dependent jet cross section at LO and NLO, 
again at y/S = 200 GeV and = 500 GeV as in Figs. 0] (a) and (b), respectively. In each case 
the shaded bands indicate the uncertainties from varying the unphysical scales in the range 
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Figure 6: As in Figs. El (a) and (b) but now showing the spin asymmetry at NLO using the GRSV 
"standard" set [25j as well as three other sets with very different gluon polarizations (see text). 
The "error bars" indicate the expected statistical accuracy beam polarization and 

integrated luminosities of 3 pb~^ and 20 pb"^ for = 200 GeV and VS = 500 GeV, respectively 
(see text). 

< /iR = /iF < 2p^T*- The solid lines are for the choice where both scales are set to p^^*. 
Clearly, the scale dependence indeed becomes much smaller at NLO, a result that was already 
noted in [9\. We emphasize that the scale ambiguity of the single-inclusive jet cross section is 
somewhat smaller than for the corresponding single- inclusive hadron cross section, cf.. Fig. 5 in 
This is not entirely unexpected, as the additional final-state factorization scale fip related 
to the fragmentation of a parton into the observed hadron provides a further source of scale 
dependence not present for jet production. 

Next we consider the double spin asymmetry v4ll' -^1- (I22I); which is the main quantity of 
interest in experiment. Fig. IHl shows A^L' calculated at NLO for the "standard" set of GRSV 
parton densities and for three other sets emerging from the GRSV analysis, which mainly 
differ in the assumptions about the gluon polarization: "A(7 = g input" , "A^f = input" , and 
"A^f = —g input" . These are characterized by a large positive, a vanishing, and a large negative 
gluon polarization, respectively, at the input scale of the GRSV analysis [21]. We should note 
that all sets provide a good description of all presently available data on spin-dependent deep- 
inelastic scattering. Again we show results for both c.m.s. energies relevant for RHIC; the 
other parameters are chosen as before. Also shown is the expected statistical accuracy for such 
measurements in certain bins of p!^^ for the Star experiment, calculated from P 

SAf' ^ } . (34) 

We assume a moderate proton beam polarization Vp of 40%. For the c.m.s. energy = 
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Figure 7: Relative contributions from gg, qg, and qq scatterings to the NLO polarized cross section 
for the "standard" set of GRSV (2H] for jet and inclusive-7r*^ production at mid pseudo-rapidities. 
Note that we use two separate coordinate axes for the two cases. 



200 GeV we use an integrated luminosity C of only 3 pb~ which should be well accomplishable 
in the near-term future. Also, we assume that only half of the calorimeter is instrumented, i.e., 
we integrate only over < rj^^^ < 1. The error bars at = 500 GeV refer to £ = 20pb~^ 
and —1 < r^j'^* < 1. Statistical accuracies for higher beam polarizations and/or luminosities are 
easily obtained by a proper rescaling of the error bars according to Eq. (j34j) . 

First and foremost we conclude from Figs. El that there are excellent overall prospects for 
determining Ag from measurements at RHIC: the spin asymmetries for the different sets 
of polarized parton densities, which mainly differ in the gluon density, show marked differences, 
much larger than the expected statistical errors in the experiment even for the very moderate 
luminosities and beam polarizations assumed here. This finding was already highlighted in [H]. 
It is interesting, however, that at moderately small p^^* the spin asymmetry is insensitive to 
the sign of Ag. A large negative gluon polarization yields an asymmetry somewhere in between 
the ones obtained with moderately and large positive gluon polarizations. This observation is 
very similar to the one recently made for inclusive hadron production at small transverse mo- 
mentum . It is due to the fact that at moderate pt and mid pseudo-rapidities the gluon-gluon 
initiated subprocess is dominant. This process has a positive analyzing power and, since the 
two gluons are predominantly probed at very similar momentum fractions, it essentially probes 
the square of Ag{xa — Xb). Thus a positive asymmetry is obtained even for a negative Ag. 
At larger jet transverse momenta, the qg process gradually takes over, resulting in sensitivity 
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Figure 8: Cone radius R dependence of the polarized and unpolarized single-inclusive jet cross 
section and the corresponding spin asymmetry at NLO for = 200 GeV and three different pj.^ 
bins. We have integrated over —1 < 7]^^^ < 1. 

of All ^'^ ^^S^ '^^ ^g{xa — Xb). When the c.m.s. energy is higher, the onset of the qg 
dominance occurs at higher p-l^^ since the spin asymmetry roughly scales with x-J^* = 2p!^^/ ^/S. 

These properties of the partonic scatterings are exemplified by the solid lines in Fig.[7[ which 
show the relative contributions to the polarized NLO cross section for the GRSV "standard" 
set for gg, qg, and qq scatterings at = 200 GeV. Here a "g" stands generically for the 
appropriate sum of all quark and anti-quark flavors. Each curve has been normalized to the 
full NLO cross section, so that the three lines add up to unity at every pl^. For illustration, 
we also show in Fig. [3 the corresponding results for inclusive 7r° production E] at mid pseudo- 
rapidities. We see that all curves for jets and pions are almost congruent, provided we rescale 
the axis for pJ^ by about a factor of 2. This feature is understood from the fact that pions result 
from a fragmentation process, in which the pion inherits only a certain momentum fraction z of 
a final-state parton. At RHIC energies, mid pseudo-rapidities, and for the transverse momenta 
we are considering here, one finds that the average z is about 0.5. This means that a pion of, say 

= 5 GeV on average originates from a scattering in which a 10 GeV parton was produced. 
For the jet cross section, this parton would produce however a jet with p!^^ = 10 GeV. This 
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explains the results in Fig. [7| Similar relations are also found for the spin asymmetries ^i-nd 
All. This interesting predicted interplay between hadron and jet observables may be exploited 
to cross-check results and to gain a deeper insight into the dynamics producing high-transverse 
momentum final states. 

We close by showing in Fig. |H1 the dependence of the polarized and unpolarized single- 
inclusive jet cross sections and the corresponding spin asymmetry Aj^L the cone size R 
in NLO QCD, for three representative bins in the jet's transverse momentum p-J^*. We have 
chosen y/S = 200 GeV and again integrated over pseudo-rapidity, — 1 < z/-'^* < 1. The results 
for = 500 GeV are very similar and not shown here. We recall from Sec. II that for R not 
too large, the dependence on R is logarithmic. Obviously, at least in the unpolarized case, the 
cross section has to rise with increasing R. At larger R, deviations from our curves would be 
expected due to the terms oc R^ becoming important. A striking finding in Fig. IHlis that All 
turns out to depend only very mildly on R, in particular for pl^^ not too large. We also note 
that since the cross section in Born approximation is independent of R, a measurement of the 
cone size dependence is a direct probe of NLO contributions. 



IV Summary and Conclusions 



In this paper we have presented a NLO calculation for the spin-dependent hadroproduction 
of single-inclusive jets. The application of the small cone approximation not only allowed 
us to perform the calculation at a largely analytical level but also to use major parts of a 
previous calculation for single-inclusive hadron production. We have outlined in some detail 
the connection between these two cross sections. By comparing to the Monte-Carlo jet code of 
[Hj which treats the cone size exactly, we have demonstrated the applicability of the SCA up to 
cone sizes of about 0.7. Our code has the advantage of being numerically very stable and fast. 

Our results are useful for studies of large-pT jet production at RHIC, in particular for 
the analysis of upcoming data in terms of the unknown spin-dependent gluon density in the 
nucleon. We found that the NLO corrections to the (un)polarized cross section are well under 
control, and that their inclusion leads to a significant reduction of scale dependence, the main 
source of theoretical ambiguity. Even for rather moderate beam polarizations and integrated 
luminosities, the corresponding spin asymmetry shows a sensitivity to A^f. This makes single- 
inclusive jet production an excellent tool for fulfilling the short-term goal of RHIC: a first 
determination of the gluon polarization. Also in the long-run jet production data will provide 
invaluable information in a detailed mapping of the x-shape of spin-dependent parton densities. 
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The experimentally relevant spin asymmetry turned out to be rather insensitive to the precise 
definition of the jet as well as to the actual choice of the cone opening. 
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